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Cometary water expansion velocity from OH line shapes 
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ABSTRACT 



Aims. We retrieve the H2O expansion velocity in a number of comets, using the 18-cm line shapes of the OH radical observed 
with the Nancay radio telescope. 

Methods. The H2O velocity is derived from the large base of a trapezium fitted to the observed spectra. This method, which 
was previously applied to 9 comets, is now extended to 30 further comets. This allows us to study the evolution of their water 
molecule outflow velocity over a large range of heliocentric distances and gas production rates. 

Results. Our analysis confirms and extends previous analyses. The retrieved expansion velocities increases with increasing gas 
production rates and decreasing heliocentric distances. Heuristic laws are proposed, which could be used for the interpretation 
of observations of cometary molecules and as a touchstone for hydrodynamical models. The expansion velocities retrieved 
from 18 cm line shapes are larger than those obtained from millimetric observations of parent molecules with smaller fields of 
view, which demonstrates the acceleration of the gas with cometocentric distance. Our results are in reasonable quantitative 
agreement with current hydrodynamical models of cometary atmospheres. 
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1. Introduction 

Due to the small gravity field of cometary nuclei, cometary 
atmospheres are not captive and expand freely in the in- 
terplanetary medium. The expansion velocity of the atmo- 
sphere is a parameter of crucial importance for the inter- 
pretation of cometary observations and the modelling of 
cometary phenomena. This velocity is not constant; it is 
governed by the sublimation mechanism, photolytic heat- 
ing of the coma and collisions, so that case-by-case studies 
are in principle necessary. It is suspected to depend basi- 
cally upon the heliocentric distance, which governs pho- 
tolytic heating, and the gas production rate, which gov- 



rectly trace the velocity distribution over the line of sight 
This has been applied to the millimetric and submillimet- 
ric lines of molecules, such as HCN, CO, CH3OH, H2S 
direc t ly sublimated from nucleus ices (e.g., |Despois et al 



1986; ISchloerb et alJll987t iBiver et alj |2002) 



erns collisions (e.g ., IBockelee-Morvan fc Crovisier! 11987 



ICombi et al. 2005 ). Heuristic laws have been proposed, 
but they need to be checked and validated by hydrody- 
namical models and direct measurements. 

One of the best methods to determine this expansion 
velocity is the observation of the shapes of molecular lines 
at radio wavelength, taking benefit of the high spectral 
resolution of such observation. These lines, when optically 
thin (which is practically the case for all species except 
water), have purely Doppler profiles. Therefore, they di- 



One can also use the OH lines at 18 cm, for which an 
important cometary database now exists. However, OH is 
a secondary species, coming from the photodissociation of 
water. At photodissociation, the O H radical is giv en an 
ejection velocity Vd w 1 km s _1 ( Crovisier 1989f ). The 
kinematics and the space distribution of the OH radi- 
cal must be evaluated in t he frame of a vectorial model 
(|Combi fc Delsemmel[l980l;lFestoulll98ll) . 

First attempts to deconvolve the OH line 
profiles in order to retrieve the expansion ve- 

Carlo 



1984 
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locity V p of t he OH-parent used Mo nte 

simulations ( Bockelee-Morvan fc Gerardl 

iTacconi-Garman et al.l Il99dh . Another approach, pro- 
posed by IBockelee-Morvan et al. (|1990L hereafter Paper 
I), derived the H2O velocity from the large base of a 
trapezium fitted to the observed 18-cm line shapes of 
the OH radical. This method, which is much more rapid 
than Monte Carlo simulations, was applied to 9 comets 
observed with the Nangay radio telescope and checked 
over the Monte Carlo method in specific cases. 
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Fig. 1. A sample of OH spectra with their fitted trapezia: 
C/200f A2 (LINEAR). 

In the present paper, we extend this work to 30 further 
comets recently observed at Nangay, including C/1995 01 
(Hale-Bopp) which was observed up to more than 4 AU 
from the Sun, and several comets which were observed 
close (s=s 0.5 AU) to the Sun. This allows us to study the 
evolution of their water molecule outflow velocity over a 
wide range of heliocentric distances and gas production 
rates. 

Section [2] presents the database and the method of 
analysis. Section [3] presents the results and the correla- 
tion between the gas expansion velocity, the heliocentric 
distance and the gas production rate. In Section HJ these 
results are discussed in the frame of hydrodynamical mod- 
els of cometary atmospheres and are compared with other 
determinations of cometary expansion velocities. Section[5] 
conclu des. Further details on this work are given in lTsengl 
(|2004h . 



2. Analysis 

The database of OH obser vations of come t s at N angay 
up to 1999 is described in Crovisier et al. (2002a '^1. In 
1995-2000, the Nanga y radio telescope was upgraded 
( van Driel et al. 19961 ). resulting in improved perfor- 
mances of the instrumentation (the sensitivity was im- 
proved by a factor « 2 and a more versatile spectrometer 
was installed). C ometary observations a fter this upgrade 
are described bv I Crovisier et al ] (|2002bl . and in prepara- 
tion). At 18 cm wavelength, the Nangay radio telescope 
has an elliptical field of view of 3.5' x 19', which corre- 
sponds to n x r 2 = 76 000 x 415 000 km at A = 1 AU 
and is equivalent to a circular field of view of radius 
r = Vn^2 = 178 000 km. 



Table 1. Comets considered for the present analysis, 
listed by order of perihelion dates. 



Comet 


N 


ret. 


C/1982 Ml Austin 


2 


a) 


C/1984 Nl Austin 


2 


a) 


21P/Giacobini-Zinner (1985) 


5 


a) 


C/1985 Rl Hartley-Good 


4 


a) 


C/1985 Tl Thiele 


1 


a) 


IP/1982 Ul Halley 


26 


a) 


C/1986 VI Sorrells 


1 


a) 


G/ 1986 PI Wilson 


13 


a) 


C/1987 PI Bradfield 


3 


a) 


C/1990 Kl Levy 


15 


b) 


109P/1992 S2 Swift-Tuttle 


10 


c) 


C/1988 Al Liller 


1 


d) 


23P/1989 Nl Brorsen-Metcalf 


2 


d) 


C/1989 Ql Okazaki-Levy-Rudenko 


2 


d) 


C/1989 Wl Aarseth-Brewington 


1 


d) 


C/1989 XI Austin 


3 


d) 


C/1991 Yl Zanotta-Brewington 


1 


d) 


C/1991 T2 Shoemaker-Levy 


1 


d) 


24P/Schaumasse (1993) 


1 


d) 


C/1993 Yl McNaught-Russell 


1 


d) 


19P/Borrelly (1994, 2001) 


8 


d) 


45P/Honda-Mrkos-Pajdusakova (1996) 


1 


d) 


C/1996 B2 Hyakutake 


9 


d) 


22P/Kopff (1996) 


2 


d) 


C/1996 Ql Tabur 


1 


d) 


C/1995 01 Hale-Bopp 


28 


d) 


C/1998 Jl SOHO 


1 


d) 


/ VjiacoDini-Zjinner (^±yyo j 


O 
Z 


A\ 
a ) 


C/1999 HI Lee 


7 


d) 


C/1999 N2 Lynn 


1 


d) 


C/1999 Tl McNaught-Hartley 


3 


d) 


C/2000 Wl Utsunomiya- Jones 


1 


d) 


C/2001 A2 LINEAR 


9 


d) 


C/2000 WMi LINEAR 


3 


d) 


153P/2002 CI Ikeya-Zhang 


8 


d) 


C/2002 Fl Utsunomiya 


2 


d) 


C/2002 VI NEAT 


6 


d) 


C/2002 X5 Kudo-Fujikawa 


3 


d) 


C/2002 Yl Juels-Holvorcem 


3 


d) 



See also 

http : //www. lesia. obspm.fr/planeto/cometes/basecom/ 



N is the number of sampl e s for each comet. 

a) iBockclcc-Morva n et alj (11990). 

b) iBockelee-Morvan et alj (Il992l ). 

c) IBockelee-Morvan et alj (j 1994T) . 

d) Present work. 



As explained in Paper I, a symmetric trapezium is fit- 
ted to the OH line. The large base of the trapezium is 
assumed to be 2(V P + Vd), where V p is the water expan- 
sion velocity and Vd is the OH ejection velocity upon wa- 
ter photodissociation. As discussed in Paper I, we assume 
Vd = 0.9 km s~ \ which is clos e to the theoretical value 
of 1.05 km s _1 (Crovisier 1989). Hence the evaluation of 
V p . An example of OH spectra with their fitted trapezia 
is shown in Fig. [1] 
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Fig. 3. Correlation between the retrieved water expansion velocity V p , heliocentric distance r^, and OH production rate 
Qoh- The size of each circle is proportional to V p . The well-studied comets lP/Halley, C/1990 Kl (Levy), 109P/Swift- 
Tuttle, C/1996 B2 (Hyakutake), C/1995 01 (Hale-Bopp) and 153P/2002 CI (Ikeya-Zhang) are identified by their 
names. 



0.5 - 



C/1995 01 Hale-Bopp 

OH observations 
Power-law fit from OH 
Power-law fit from millimetric data 



f--~ - I:? - J 



1 2 3 4 5 

Heliocentric distance [AU] 

Fig. 2. Evolution of the expansion velocity V p as a func- 
tion of heliocentric distance r/,. for C/1995 01 (Hale- 
Bopp). Pre-pcrihclion and post-perihelion data are com- 
bined. The dashed line shows the power-law fit (Eq. [IJ 
to the OH data. The dotted line shows the power-law 
fit (Eq. [J) to t he velocities deriv ed from the millimetric 
molecular lines (jBiver et al.ll2002l ). 



As in Paper I, the OH producti on rates Qqh are de - 
rived using the OH inversion curve of lDespois et al I (Il98lh 



and t he quenching law and OH param eters of | Ge rard 



1990). The model is fully described in ICrovisier et al 



2002a, Table 3, last column). The model consistently 



uses the parent velocity derived from the trapezium fit 
to compute the OH production rate; therefore, the pro- 
duction rates of the prese nt analysis may d i ffer som ewhat 
from those published by ICrovisier et al. (|2002al ) which 
were computed assuming V p — 0.8 km s . 

The list of comets investigated in the present study 
is given in Table [1] The analysis is performed on spectra 
integrated over several days for which the signal-to-noise 
ratio is sufficient (typically > 10). The total number of 
samples is 190. A comprehen sive tabulation of the data 
may be found in iTsen K (2004). 

As an example, the expansion velocity as a function of 
heliocentric distance is shown in Fig. [2] for C/1995 01 
(Hale-Bopp), for which the OH lines could be observed 
over a large range of r^. The fitted power law is 



V p = 1.917(±0.055> 
to be compared to 



-0.78(±0.04) 



km s 



(1) 



V p = 1.125(±0.015)r~ 



0.42(±0.01) 



km s 



(2) 
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1 .0x10 < OH production rate < 1 .0x10 s 



0.3 AU < heliocentric distance < 0.6 AU 



2.5 - 



0.5 




12 3 4 

Heliocentric distance [AU] 

1 .0x1 29 < OH production rate < 1 .0x1 30 s" 1 



0.5 




2 3 
Heliocentric distance [AU] 



lines (jBiver et al.ll2002i ) which were observed in a smaller 
field of view (see discussion in Section d]). One must note 
that in the range 0.92 < < 4.57 AU which is investi- 
gated here, Qcm varied by about three orders of magni- 
tude, so that the strong variation on implicitly includes 
a dependence on both rh and Qoh- 

3. Results 

From the profiles of the 18-cm OH lines, we have investi- 
gated the expansion velocity V p of cometary atmospheres 
for heliocentric distances r/j ranging from 0.3 to 4.6 AU 
and OH production rates Qoh from 10 28 to 10 31 s _1 . 
Confirming the results of Paper I, V p is found to consis- 
tently increase with increasing Qoh and with decreasing 
Th„ as expected qua litatively from hydrodynamical models 
( Combi et al. 120051 and references therein). The observed 
Vp's range from 0.5 to 2.5 km s _1 . 

An overview of the correlation between retrieved water 
expansion velocities, heliocentric distances and OH pro- 
duction rates is shown in Fig. [3] (an update of Fig. 14 in 
Paper I). In order to investigate separately the dependence 



2.5 



1.5 



0.5 



Fig. 4. Evolution of the expansion velocity V p as a func- 
tion of heliocentric distance rh for given Qoh ranges. The 
fitted power laws are listed in Table O 



for the velocities derived from the millimctric molecular r- 
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Fig. 5. Evolution of the expansion velocity V p as a func- 
tion of Qoh for given heliocentric distance (?h) ranges. 
The fitted power laws are listed in Tabled 



of V p to r/,, and Qoh, we have divided the data into various 
subsamples. Power laws fitted to the different subsamples 
are listed in Table and shown in Figs 2H5] 

From these plots, we can see the existence of a thresh- 
old effect: the expansion velocities are « 0.8 km s _1 , in- 
sensitive to the heliocentric distances from 1.5 to 5.0 AU 
for moderately active comets with OH production rates 
from a few 10 28 s _1 to a few 10 29 s _1 . The expansion ve- 
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Fig. 5. (continued) 



locities strongly depend on the heliocentric distance for 
r h < 1.5 AU. 

Errors reported in Table [2] are the 1-tr formal errors 
resulting from the fits to the data. They do not include 

Table 2. Power-law fits to the H2O expansion velocity V p 
(in km s~ x ) obtained from the OH radio observations at 
Nancay for different ranges of production rates Qoh and 
heliocentric distances r^. The corresponding data with 
plots of the fitted power laws are shown in Figs 0H5] 



Qoa range 
Is" 1 ] 




power-law fit 
[rh in units of AU] 


10 28 _ 10 29 


v P 


= 0.96 ±0.01 x r h -°- 44±0 ' 04 


10 29 - in 30 


v P 


= 1.19 ± 0.01 x r ;i - - 55±0 ' 02 


r h range 
[AU] 




power-law fit 
[Qoh in units of 10 29 s" 1 ] 


0.3 - 0.6 


Vp 


= 1.42 ± 0.05 x <2oH a30±0 ' 06 


0.6 - 0.9 


Vp 


= 1.23 ± 0.02 x Qoh°' 19±0 ' 01 


0.9 - 1.1 


Vp 


= 1.11 ± 0.02 x Qoh°- 22±0 ' 01 


1.1 - 1.3 


Vp 


= 1.07 ±0.02 x Qoh 015±0 ' 01 


1.3 - 1.5 


Vp 


= 0.78 ± 0.02 x QoH 019±0 ' 02 



possible systematic errors due to shortcomings in the mod- 
elling. Given the limited number of data, their arbitrary 
binning, and possible systematic effects, we do not think 
that the differences of the exponent of the power law r%, 
in the upper part of Table [2l are firmly established. The 
same for the exponent of Q Q h m the l° wer part of the 
table, a = —0.5 and /3 = 0.2 are probably the values to 
be retained. 

Taking a step further, we tried to fit a combined power 
law r^Qo H t° our data. A unique law could not be ob- 
tained for the whole domain covered by the observations. 
The results are shown in Table [3] For rh > 2 AU, V p barely 
depends upon rh and Qoh as already discussed. The de- 
pendence with Qoh is oc Qoh° over * ue ra nge 0.3-2 AU, 
as inferred from Qoh power-law fits over restricted ranges 
of rh (Table[2|). Between 0.3 and 1 AU, the dependence 
is in r^ ' 45 , but the data suggest a steeper heliocentric de- 
pendence in r^ 1 between 1 and 2 AU. Power-law fits in the 
subsamples [10 28 < Q OH < 10 29 s" 1 ; 1 < r h < 1.5 AU] 
and [10 29 < Qoh < 10 31 s" 1 ; 1 < r h < 2 AU] are similar 



oh ' h 



which suggests that the difference in helio- 
centric evolution between < 1 AU and rh > 1 AU could 
be a characteristics of V p behaviour. Note, however, that 
there is significant dispersion between the data and the 
fitted curves (reduced x 2 between 2 and 3, see Table 



4. Discussion 

The empirical laws derived from non-linear fits and listed 
in Table[5]can be used as a touchstone for hydrodynamical 
models of cometary atmospheres. They can also help us 
to analyse cometary observations, and especially to derive 
more reliable OH production rates from the 18-cm obser- 
vations, when the signal-to-noise ratio is not sufficient to 
apply the trapezium method to the line shapes. 

These results, however, should only be used for the 
range of parameters (rh, Qoh) pertaining to the data 
analysed in the present work. In particular, we caution 
the reader against extrapolating the laws of Table [2] to 
distant comets. At large r/i's, cometary activity is gov- 
erned by the sublimation of CO and other hypervolatile 
species rather than that of water, resulting in quite differ- 
ent t emperatures and velocities at the nucleus surface flpl 
1983h . 

The top part of Table [5] shows that, for a restricted 
range of Qoh, the heliocentric dependence is close to the 
law 



V v = 0.85 x r 



-0.5 



km s 



(3) 



which was proposed some time (e.g., Cochran fc Schleicherl 



1993; Budzicn ct a 



cannot be ignored, 
the law 

V v = 0.7+1.3 x 



lll994 ). But the dependence upon Qoh 



Cochran k. Schleicher ( 19931 ) proposed 



?h 2 o°' 5 km s~ 



(4) 



for rh ~ 1 AU from the results of Paper I, where Qh 2 is 



in units of 10 30 
our data. 



This law appears to be a poor fit to 
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Table 3. Combined 77, and Qoh power-law fits to the expansion velocity V p (in km s _1 ) for different ranges of Qoh 
and heliocentric distance. 

ru range Qoh range combined power-law fit reduced \ 2 

[AU] [s- 1 ] [Qoh in units of 10 29 s" 1 ] 

0.3 < r h < 1.0 1 28 < Qoh < 10 30 V v = 1.11 ± 0.02 x r -°-«±o.o5 x Q o.23±o.oi ~ 

1.0 <r h < 2.0 10 28 < Qoh < 10 31 V p = 1.17 ± 0.02 x r -°- 99±0 - 0B x Q o ^ 8±0 - 01 3.3 

2.0 <r h < 4.6 1 28 < Qoh < 10 30 V p = 0.71 ± 0.02 x r -° 07 ± 08 x Q o ° 9±0 02 2.0 



Interestingly, photochemical heating is more sensitive 
to the heliocentric distance than to the water produc- 
tion rate. The heating rate is directly proportional to 
r^ 2 (through the water photodissociation rate), but is 
weighted by the efficiency of the thermalization process 
of fast hydrogens. The size of the region where the effi- 
ciency is significant sc ales proportionally to the water pro 



ducti on rate (e.g., 13 1983 : Bockelee-Morvan fc Crovisier 



19871 ). Though other processes, such as radiative cooling, 
do affect the hydrodynamics of the coma, the steeper rh 
dependency (compared to Qoh variation) observed in our 
data is consistent with photolytic heating being the main 
process controlling the gas velocity in the outer coma. The 
different heliocentric variations in the 0.3-1 AU and 1-2 
AU domains remain to be explained. 

We will now compare our results with other observa- 
tions and with model predictions. For this prospect, we 
have to take into account the effect of the field of view 
and of the geocentric distance A. 

The expansion velocity of a cometary atmosphere is 
not a well defined parameter. In the inner, collisional coma 
where classical hydrodynamics prevails, this velocity is 
progressively increasing with dist ance r to the nuc leus, 
as a result of photolytic heating ( Combi et al.l [2005^ . In 
the outer coma where the free-molecular flow is governed 
by rarefied gas dynamics, the velocity distribution is not 
maxwellian. 

The observed line shapes are averages of all molecules 
present in the instrumental field of view. In the following, 
we will assume that the observed V p is representative of 
molecules at a distance r equal to the field-of-view radius. 
This assumption is justified by Monte Carlo simulations of 
the HCN line shape (Paper I, Fig. 15). Alternatively, the 
field of view radius for 18-cm observations may be larger 
than the water scale length l p for photodissociation (typ- 
ically l p ss 10 5 km at 1 AU). In this case, we may wish to 
adopt r = l p . This must be considered as an approxima- 
tion. 

We note that one would also expect that for very small 
fields of view, the collisional region is sampled, in which 
OH is thermalized and partakes the kinetics of water. 
Then the trapezium model would fail. An extreme obser- 
vational case was that of comet C/1996 B2 (Hyakutake), 
which was observed at = 1.06 AU, A = 0.13 AU, with 
Qoh w 3 x 10 29 s _1 . We retrieved V p = 1.4 km s _1 , right 
in the range of the laws of Table [D Quite recently, we ob- 
served 73P/Schwassmann-Wachmann 3 at A = 0.08 AU 



only, with r h = 1.0 AU and Qoh ~ 10 28 s _1 ; the re- 
trieved V p (0.8 km s" 1 ) was also normal (Crovisier et al., 
in preparation; these observations are not in the database 
analysed in the present work) . This suggests that for both 
observations, the field of view was still not small enough 
for a significant thermalization of the OH radicals, or that 
this thermalization also quenched the OH maser, so that 
the 18 cm observations are insensitive to such OH radicals. 

Complementary data are provided by millimetric 
and submillimctric observations of line shapes of par- 
ent molecules. Vp is then close to the line half-width 
at half maximum ( Biver et al. 20021 adopt for V p 90% of 
the half- width, to account for thermal broadening). Data 
for several molecules (especially HCN, CO, CH3OH and 
H2S) are now available for many com ets (e.g., observa- 
tions at IRAM, S EST, JCMT and CSO: lBiver et al.lll999l . 
2000l 120021 l2006h . Direct observations of the water line at 



557 GHz from space with the SWA S and Odin sate llites 
also exist (jNeufeld et al.ll200"fl Lecacheux et al.ll2003 ). but 
their kinematic interpretation is hampered by the satura- 
tion of this strong line. 

The typical angular field-of-vicw diameter for millimet- 
ric observations with the IRAM 30-m telescope is 17" at 
a frequency of 145 GHz, corresponding to a linear radius 
r = 6 200 km at a geocentric distance A = 1 AU. This is 
much smaller than the elliptical field of view of the Nangay 
radio telescope (equivalent to a radius of 178 000 km). 

As previously noted in, e.g., Paper I, the V p determined 
from parent molecules observed at millimetric wavelengths 
are systematically smaller than those obtained from OH 
line shapes, with larger fields of view. This effect is shown 
in Fig. IB] Observations of comets lP/Halley, C/1996 B2 
(Hyakutake), C/1995 01 (Hale-Bopp), C/1999 HI (Lee) 
and 153P/2002 CI (Ikeya-Zhang) were selected at vari- 
ous tVs, chosen so that the gas production rates were 
similar (Qoh in the range 1-4 xlO 29 s" 1 ). Points on the 
right are V p determined from the present work. Points on 
the left are measurements from m illimetric observations 
(|Biver et alll999l . l2000L[20oll2006h . For comet Halley, the 
various points f rom left to right are from in situ measure- 
ments (Giotto, iLammerzahl et all Il987l). HCN (IRAM . 
Despois et aDl986h . HCN (FCRAO JSchloerb et alll987h 
and OH (Nangay, present work). As discussed above, r 
is evaluated as the field-of-view radius. For all comets, 
an increase of V p with r is observed. This increase is 
much more important for comets close to the Sun (e.g., 
153P/2002 CI (Ikeya-Zhang) at r h = 0.6 AU) than for 
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1 .0x10 < OH production rate < 4.0x10 s 



2.5 



1.5 



Ikeya-Zhang (0.6 AU) 
C/1999 H1 Lee (0.85 AU) 
Halley (0.9 AU) 
Hyakutake (1 .5A U) 
Hale-Bopp (3.0 AU) 



3 3.5 4 4.5 5 5.5 6 

log nucleocentric distance [km] 

Fig. 6. Ficld-of-view effects demonstrating the accelera- 
tion of the gas. V p determined from the present work 
(right-hand points) or from millimetric observations (left- 
hand points) is plotted as a function of the cometocentric 
distance r corresponding to the field of view. See text for 
details. 



3.0 AU, all pre-perihelion) and the corresponding observed 
water production rates (80, 20, 14, 4 and 4 x 10 29 s _1 , re- 
spectively, according to Combi et al.). Observational de- 
terminations of V p , from the present work (p oints on the 
right) and from the millimetric observations of lBiver et alj 
(2002]) (points on the left), have been superimposed. For 
all these points, the cometocentric distances are again 
those corresponding to the field-of-view radii. We can see 
that the model predictions lead to somewhat overesti- 
mated velocities. Otherwise, the variation of V p with r 
and rh (and the related Qoh) is remarkably well repro- 
duced by the model. The discrepancy between model and 
observations may be attributed to: 

flaws in the modelling of the heati ng/cooling pr o cesses 



(one can note that the model of Combi et al. ( 19991 ) 



als o predicts temper atures higher than those observed 
bv lBiver et aT1l2002lk 

the difficulty to model the transitional region between 
hydrodynamical flow and free-molecular flow (s ophisti- 
cated Monte Carlo simulations are required; cf . iHodgesl 
1993); 



C/1995 01 Hale-Bopp 

1.0 AU 
1.5 AU 
2.0 AU 
2.5 AU 
3.0 AU 

millimetric data (field of view) + 
Nangay observations (field of view) x 
Nancay observations (scalelength) * 




3.5 
3 

2.5 
2 

1.5 
1 

0.5 



1 2 3 4 5 6 

log nucleocentric distance [km] 

Fig. 7. Model predictions of V p as a function of cometo- 
centric distance, for comet C/1995 01 (Hale-Bopp) pre - 
perihelion, adapted from Fig. 1 of ICombi et all (jl999h . 
Observational determinations have been superimposed. 
See text for details. 



distant comets (e.g., C/1995 01 (Hale-Bopp) at 3 AU). 
This clearly demonstrates the acceleration of the gas with 
cometocentric distance, attributed to photolytic heating. 

Direct imag i ng of near-UV OH emission was made 
bv lHarris et ail (|2002h in comet C/1995 01 (Hale-Bopp). 



— an inaccurate estimation of the cometocentric distance 
to which the measured V p pertains (see discussion 
above) ; 

— the presence of anisotropic outgassing, whereas models 
assume spherical symmetry. 

Finally, we note that, as already discussed in Paper I, 
there are some observations, especially for distant, weak 
comets, leading to small values of V p (0.5-0.6 km s" 1 ) that 
can difficultly be reconciled with hydrodynamical models. 

5. Conclusion 

We have extended the previous analysis of OH 18-cm line 
shapes from the Nangay database (paper I) from 9 to 39 
comets with the following results: 

— This analysis confirms the increase of the coma expan- 
sion velocity V p with increasing gas production rate 
Qoh and decreasing heliocentric distance r/j. 

— The results are summarized in Tables [2] and [3] which 
give power-law fits to the data. These suggested laws 
can be used for testing hydrodynamical models of 
cometary atmospheres and for the interpretation of 
molecular observations, especially OH observations 
made with similar fields of views. 

— The present results, compared to analyses of line 
shapes of parent molecules observed at millimetric 
wavelengths, yield significantly larger V p s. This may 
be attributed to the much larger field of view of the 18- 
cm observations. The two sets of data are complemen- 
tary. The increase of V p with the field of view demon- 
strates the acceleration of cometary gas with cometo- 
centric distance. 

Bopp) . V p as a function of cometocentric distance is plot- — A threshold effect is observed: comets with small gas 
ted for five heliocentric distances (1.0, 1.5, 2.0, 2.5 and production rates (typically Qoh < a few 10 28 s _1 ) 



From their analysis, V p — 2.3 km s _1 was retrieved for 
Qoh ~ 10 31 s™ 1 for cometocentric distances « 10 6 km. 
This agrees well with our measurement (2.4 km s _1 ). 

In order to achieve a quantitative comparison of these 
results with models, we have plotted in Fig. [7] the pre- 
dictions of the hybrid kinet ic /dusty gas hydrodynamical 
model of lCombi et all (|l999h for comet C/1995 01 (Hale- 
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have similar expansion velocities (typically V p « 0.8 
km s" 1 ). Gas acceleration is inefficient for such weak 
comets. 

— Our results are in reasonable agreement with cur- 
rent hydrodynamical models of cometary atmospheres. 
Small discrepancies may be attributed to modelling 
issues such as inadequate treatment of the cool- 
ing/heating processes, or of the transitional region be- 
tween hydrodynamical flow and free-molecular flow. 
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